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Abstract 

We carry a Monte Carlo study of the coupled two-scalar Xcfyfcf)^ 
model in three dimensions. We find no trace of Inverse Symmetry 
Breaking in the region of negative A's for which the one-loop effec- 
tive potential predicts this phenomenon. Moreover, for A's negative 
enough, but still in the stability region for the potential, one of the 
transitions turns out to be of first order, both for zero and finite tem- 
perature. 
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1 Introduction 



It is well known today that the strong and electroweak interactions at the low 
temperatures and energies characteristic of the present-day Universe are well 
described by a spontaneously broken relativistic gauge theory. According 
to this picture, a region of empty space resembles a ferromagnet below the 
Curie temperature, since it is characterized by a number of non-vanishing 
order parameters (one or more, depending on the elementary particle model) 
which break some of its internal symmetries, an analogue of spontaneous 
magnetization for ferromagnets. Another piece of wisdom, from statistical 
mechanics this time, is that an ordered system, when its temperature is 
raised, gradually looses its order undergoing one or more phase transitions 
until it reaches a temperature above which no order at all is present anymore. 
Long ago Kirzhnits and Linde argued that something similar might have 
occurred in the early Universe, when high temperatures were present, so that 
the internal symmetries that appear broken today might in fact have been 
manifest at that time. Since then this intuitive picture has been confirmed 
by detailed computations [TJ, |3| and computer simulations. The idea that 
the ground state of the Universe was symmetric in the early times and that 
later, when the Universe cooled down during its expansion, there took place a 
series of phase transitions that eventually led to the present day asymmetric 
vacuum, is of the greatest importance for the structure of the Universe as it 
influences such important issues as baryogenesis, the formation and dynamics 
of topological defects [f| like monopoles and cosmic strings, just to mention 
a few. 

Is this scenario the only possible one? Maybe not. In fact, already in the 
classic paper by S. Weinberg |2[ it was pointed out that the degree of order 
of the vacuum may increase when the temperature is raised, contrary to any 
intuition, in models with a sufficiently reach scalar sector. This phenomenon 
was called Symmetry Non-Restoration (SNR) or Inverse Symmetry Breaking 
(ISB), depending on whether the vacuum is ordered or disordered at zero 
temperature. It is due to the possibility that some of the scalar fields acquire 
a negative Debye mass, via negative quartic interactions among themselves 
and with the other scalars, something which is allowed in multiscalar models 
without causing any instability in the potential. To make this seemingly 
paradoxical statement more plausible, the author observed that in Nature 
there exists at least one substance, the Rochelle salt, which does exhibit 
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this type of "inverse" behavior in a certain range of temperatures and thus 
the expectation that more heat necessarily means more disorder is clearly 
not always true. In any case, these remarks of Weinberg passed totally 
unnoticed, until some authors resumed this idea, applying it to important 
cosmological questions like the domain wall and monopole problems || 
0, ||, the breaking of the CP symmetry ||, baryogenesis |I0| , inflation |fLl] 



and the breaking of the P, Strong CP and the Peccei-Quinn symmetries [fL2 
Recently, the possibility of ISB and SNR also in supersymmetric models has 
been explored Jl3| . 

Despite the potentially important consequences for Cosmology, in the 
literature on ISB and SNR there is a certain amount of scepticism about the 
very existence of these phenomena. The example of the Rochelle salt is not 
really persuading, because it is obvious that if a piece of it is heated enough 
it eventually melts, while in a system with ISB or SNR there is order at 
arbitrarily large temperatures. Moreover the result of Weinberg was based 
on a one-loop perturbative computation and so there is the possibility that 
it is an artifact of this approximation. Since then, many authors have tried 
to improve the one-loop result, by resorting to a number of techniques all 
including some amount of non-perturbative physics. It seems fair to say that 
the current situation is perplexing: while some studies have confirmed that 
ISB and SNR exist |15], [U| [IT], |T8fl , even though the region of parameter 



space for which they occur appears to be reduced in size compared to the 
lowest-order result, others have come to the conclusion that these phenomena 



disappear once the non-perturbative information is put in [19, EQj. In any 



case, no one of these results appears as conclusive, because all of them rely on 
approximations, with various degree of applicability. The only exact result 
known so far regards the lattice, where it can be shown that no order is 
possible at sufficiently high temperatures for arbitrary models, gauge or not 
Unfortunately these theorems have been proven only for lattices with 
finite spacing in the space-directions (but with continuum euclidean-time 
axis) and their impact for the continuum is thus unclear. 

In this paper we present the results of the first Monte Carlo study of ISB. 
We have considered a two-scalar model in 2+1 euclidean dimensions, for 
which a one-loop computation predicts the possibility of ISB and SNR. We 
have studied its phase diagram for various choices of the coupling constants 
as a function of the temperature. Concretely, we have simulated the model 
on various asymmetric lattices N t x N^, where N t and N s are respectively the 



3 



number of sites in the "euclidean time" and space directions, the temperature 
being related to N t by T = l/(aN t ), with a the lattice spacing. In order to 
avoid finite-size effects we have taken the thermodynamic limit N s —* oo (for 
any fixed N t ) using Finite Size Scaling (FSS) techniques. The high level of 
precision required to clearly separate the transitions corresponding to the 
various values of N t required very large statistics and correspondingly long 
simulation times. We have found that for all the values of N t that we have 
considered, ISB seems ruled out, irrespective of the values of the couplings. 
In fact, it appears that the size of the disordered region of the phase diagram 
increases when the temperature is raised (i.e. when N t is reduced), as it 
happens in normal cases (for example in the <p\ model J22|) , and does not 
decrease as required for ISB to take place. As for SNR the results of our 
simulations show that, if one starts at T = from an ordered phase, the 
system disorders above a certain critical temperature T c (coherently with the 
results of [PUD, again as it happens in normal cases. In principle, this does 
not exclude the possibility of SNR, because it might happen that T c diverges 
in the continuum limit. In order to see if this possibility really occurs, one 
would have to carry a detailed study of the continuum limit, something that 
we have not done. 

Another interesting result of our simulations is that one of the transitions 
turned out to be of first order, both at T = and at finite temperature, 
when the coupling among the two fields is negative and large in absolute 
value, being of second order otherwise. This result is in qualitative agreement 



with refs.[23, 24], where our model is considered for the particular case when 
an extra 0i <-> 02 symmetry is present (this case of course excludes the 
possibility of ISB or SNR). Using the method of the Average Effective Action 
[p6[| , the authors found in the phase diagram surfaces of first order phase 
transitions, but there this occurred for all negative values of the coupling 
among the two scalars, while we seem to find this behavior only for couplings 
strongly negative. 

The organization of the paper is the following. In Sec. 2 we present 
the continuum model, and discuss the main features of its lattice version, 
while in Sec. 3 we review the predictions of perturbation theory on its high- 
temperature behavior. The simulation, the techniques used and the results 
obtained are presented in Sec. 4. Finally Sec. 5 is devoted to the conclusions. 
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2 The model and its lattice formulation 



We consider the theory for two real scalar fields in 3 euclidean dimensions, 
described by the bare (euclidean) action: 

S = J d'x{\ E [(d^if + \mT&i + %*f\ + \9m}- (1) 

What will be essential for ISB and SNR, in the above action the quartic 
coupling g can be negative, as the condition of boundedness from below of 
the potential is satisfied if: 

9 < 0, 

9i 92 > 9 g 2 , g x > 0, g 2 > 0. (2) 

As it is well known scalar models with quartic couplings, like ([!]), in three 
dimensions are superrenormalizable. In the perturbative expansion of the 
zero-temperature Green's functions there is only a finite number of prim- 
itively divergent diagrams. Depending on the regularization scheme, UV 
divergencies are found only in two classes of diagrams: the tadpole dia- 
grams, at one loop, and the sunset diagrams (for zero external momenta), 
at two loops, both contributing to the self-energies. As a consequence, in 
([!]) the only quantities that require an infinite renormalization are the bare 
masses m^ 2 , while the fields and coupling constants can be identified with 
the renormalized ones. It is for this reason that in (|I|) we have written the 
latter without the suffix B. 

When regularized on an infinite three-dimensional cubic lattice of points 
Q with lattice spacing a, the above action is replaced by its discretized version 

Sl = E « 3 { E [E-CAjf^) 2 ^) + \rnT®U*) + %*U*)] + 

xen i=i,2 m / 4 - 

+f ^(x)^(x)}, (3) 

where is the lattice derivative operator in the direction //: 

AW^x) = <M* + /j*)-<M*) . (4) 
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We find it convenient to measure all dimensionful quantities in ([3]) in units 
of the lattice spacing; thus we define: 

4>i,L{ x ) = a 1/2 $i,L{x)i Ui = agt, r, = a 2 rn$ 2 u = ag . (5) 
In terms of the dimensionless quantities the lattice action now reads: 



xen i=i,2 a* 



+ 



+ 4 4>l, L {^l L {x)} , (6) 

The standard lattice notation is obtained with a further redefinition of the 
fields and couplings in (H) according to: 

</>*,!, W = V K i<Pi,r, ^ = — , Ti = 2 (7) 

and 

u = . (8) 

After these redefinitions we get our final form of the lattice action 



S L = Yl { H [ _ K i H <f>i,r<j)i,r+fi + Ai(0? r - l) 2 + 



%,r 

,. eZ 3 ' i=l,2 ' 



+ 



+ MlAl,r}- (9) 

For generic values of the parameters, the model has a Z 2 x Z 2 symmetry 
which can be spontaneously broken. We now briefly comment on the fixed- 
point structure of the model. Besides the Gaussian fixed point, corresponding 
to T\ = r 2 = U\ = u 2 = u = (or K\ = k 2 = 1/3 and Ai = A 2 = A = 0), 
which has a null attractive domain in the infrared, there are five more fixed 
points. They are all attractive in at least one direction and correspond to: 

1) the Heisenberg fixed point, for u\ = u 2 = 3u = u* H and t\ = r 2 = r* H < 0, 
where the symmetry of the model is enhanced to 0(2). 

2) the three Ising fixed points, for u = and Ui = u*j, = r} for some i, 
where the model splits into two independent <p\ models. 

3) the Cubic fixed point, for = u = u}/2, ri = r}, which again splits in 
two independent 0| models, after a 7r/4-rotation of the fields. 
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Due to the rich fixed point structure, our two-scalar model exhibits com- 
plicated cross-over phenomena. For a study of these aspects (for the partic- 
ular case when an extra <f>\ ^ 02 symmetry is present), we refer the reader 



to ref. 23 



We now briefly discuss the continuum limit of the lattice model (y), even 
if in our simulations we have not taken this step because the results we got 
seem to rule out the possibility of ISB a priori, even in the continuum. Taking 
the continuum limit would instead be necessary to prove that neither SNR 
occurs, a possibility that we cannot exclude in principle, even if it appears 
unlikely. 

In order to go to the continuum, one has to move the simulation point 
along (zero temperature) renormalization group trajectories of the lattice 
model, called curves of constant physics (CCP). They are parametrized by 
the lattice spacing a and are such that any observable, measured in physical 
units, approaches a definite limit when a — > 0. The CCP's that one has to 
pick in order to explore the models described by perturbation theory in the 
continuum are those that approach the Gaussian fixed point. This is rather 
evident if one recalls that in perturbation theory the bare coupling constants 
gi and g are finite and so, by (|5|), the lattice coupling constants Ui and u all 
vanish in the limit of the lattice spacing going to zero. We conclude these 
remarks by observing that, due to the superrenormalizability of the model, 
the asymptotic form of the CCP can be computed exactly, by means of a 
simple two-loop computation. For more details, in the case of the pure 0f 



model, we refer the reader to ref. [22 



3 High-T perturbation theory 

In this Section we briefly review the predictions of perturbation theory on the 
high-temperature behavior of the model (|I|). This turns out to be not an easy 
task. In fact, while the breaking of perturbation theory at high temperature 
is a generic feature of quantum field theories, in three dimensions things are 
worse due to the presence of severe infrared divergences. Let us see how 
this comes about. The issue we are interested in is the possibility that in 
our two scalar model the symmetric vacuum of the theory is unstable at 
arbitrarily high temperatures, for opportune choices of the parameters. For 
this we need compute the leading high-T contribution to the second derivative 
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of the effective potential in the origin, or equivalently to the 1PI 2-point 
functions for zero external momenta p. An easy one-loop computation (in 
the imaginary time formalism) |2S[ in our model gives, in the limit T/m ^> 1, 
the result: 

M?(T) = -r\ 2 \p = 0,T) = ^(g l + g)\og(^j , (10) 

where fi is the mass scale. Recalling now that the coupling g can be negative, 
it is easy to see that for (71/(72 sufficiently large, the range (H) of stability for 
the potential includes values of g such that, say, g 2 + g is negative (it can be 
proven easily form eqs.(^) that one cannot have g\ + g < at the same time). 
If we can trust eq. fllCf) , it thus appears that, for g 2 + g < 0, M|(T) becomes 
negative at sufficiently high temperatures, irrespective of its value at T = 0. 
This means that above a sufficiently high temperature the symmetric vacuum 
will become unstable in the direction of 02 arid thus there will be symmetry 
breaking for this field. This is the essence of the phenomena of SNR and ISB: 
in multiscalar models some Debye masses can become negative and so one 
can have spontaneous symmetry breaking at arbitrarily high temperatures. 

Can we trust this result in our model? As we pointed out earlier the 
perturbative expansion of the effective potential in three dimensions suffers 
from hard infrared divergencies that can cause its breakdown at high tem- 
peratures. A closer look at eq.fllOD raises the doubt that indeed this is the 
case. According to eq. (]T0|) , we see that for g 2 + g < 0, 02 acquires at high 
temperatures a vev (0 2 ) of order yT (upto logarithms) . Now, simple power 
counting arguments suggest that the loop expansion of the effective potential 
for 02 is reliable only if 



92T 



^ 2 + k 2 <02> : 



« 1 , (11) 



in the range of temperatures and fields of interest. For (02) 2 ~ T this 
quantity is of order 1 and thus it is clear that we cannot trust perturbation 
theory in this case. 

Such a conclusion is further reinforced if we look at the effect of the inclu- 
sion of higher order corrections. Indeed, at high T, powers of T can compen- 
sate for powers of the coupling constants, thus producing large corrections to 
the one-loop result (for an analysis of this problem in the (f)\ theory see the 
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first of ref . |[20|| ) . The standard solution to this difficulty is the resummation 
of the so-called ring-diagrams ||], which leads to a self-consistent gap equa- 
tion for the Debye masses Mf{T). In our case the gap equations obtained in 
this way are: 

M 2 (T) = m 2 + Z_Llog— + #log ' * 




Afi(T) * M 2 (T) J ' 
T T \ 

M 2 2 (T) = ml + -f- [ fl 2 log t-TTT^V + ^log r=-r , (12) 

2y ' 2 4vr y 6 M 2 (T) y Mi(T) / v 7 

where m? denote the renormalized masses. These equations always admit real 
positive solutions for M 2 (T), for sufficiently high T |2T) |, which means that 



the symmetric vacuum is stable at high temperatures. Thus, the inclusion 
of higher order effects (in fact the use of gap equation is a non-perturbative 
procedure, as it implies an infinite resummation) has reversed the one-loop 
result. Notice though that, differently from what happens in four dimensions 
H, the gap equations ( |T2"D do not make sense for Mf = and thus, if there 
is spontaneous symmetry breaking at T = 0, they cannot be used to give an 
estimate of the critical temperature T c above which symmetry is restored. 
This is another consequence of the infrared singularities alluded to before 
(for the case of a 0| model, the limits on the possibility of determining T c 
in perturbation theory are discussed in [^7| , while a measurement of T c as a 
function of the renormalized parameters by a Monte Carlo simulation is given 
in mi ). This state of affairs is obviously unsatisfactory: perturbation theory 
cannot give any information in the study of ISB and SNR in three dimensions 
and this motivated us to perform a lattice Monte Carlo simulation of this 
model to have a fully non-perturbative approach to this problem. 



4 The simulation 

As we saw in the previous section, according to a one loop computation of 
the Debye masses, at high temperatures the symmetric vacuum of our system 
should be unstable in the direction of, say, the field 4> 2 for negative values of 
g such that gi + g < 0. This means that, for such values of g, the field (f>2 
should necessarily have a non-vanishing vacuum expectation value (vev) at 
high temperature. 
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Let us see how one can analyze this phenomenon on the lattice. To study 
the effects of temperature on the system, one simulates it on asymmetric 
lattices N t x iV s 2 , N t and N s being the number of sites in the euclidean time 
and space directions respectively. The physical temperature is then related to 
N t as T = l/(aN t ), a being the lattice spacing. The T = case, corresponds 
to symmetric lattices N%. In the next step one takes the thermodynamic limit 
N s — > oo, holding N t fixed: for each N t , in this limit the phase diagram of 
the system splits into four phases, distinguished by the vev's of the fields and 
divided by surfaces of phase transitions. Finally, one takes the continuum 
limit a — > 0, which implies, as explained in section 2, moving the simulation 
point along some CCP approaching the gaussian fixed point. In order to see 
if ISB occurs, one now picks a CCP that describes a continuum theory with a 
symmetric vacuum at T = 0. Such a CCP must lie entirely in the disordered 
region of the T = phase diagram. One has ISB if this CCP crosses the 02 
transitions line, separating the symmetric phase from the ordered phase with 
(02) 7^ 0, of some finite Nt lattice. It is obvious that this can happen only 
if the 02 transition line of some lattice with finite Nt lies in the disordered 
region of the T = phase diagram. If a is the lattice spacing corresponding 
to the crossing point, one concludes that the physical temperature for the 02- 
ISB-phase transition is T c = l/(dN t ). In order for this transition to survive in 
the continuum limit, our CCP must in-fact cut all the 02 transitions lines for 
sufficiently large iV t 's and it must so happen that the corresponding critical 
temperatures approach a finite value for a — > 0. 

It is clear, from the above considerations, that the first thing to do in order 
to study ISB is to draw the phase diagrams of our model, for various values 
of Nt (including of course for symmetric lattices). With the parametrization 
of (|9|) the region of ISB for 02 corresponds to 

A + 6— A 2 <0. (13) 

Our model contains five independent parameters, too many for a complete 
analysis to be possible. Thus, we decided to fix the values of Ai and A 2 . 
To choose them in an optimal way, it was necessary to make a compromise 
between a number of conditions. First, we took them such that 9A2 < Ai. 
This is because the stability condition eq.(^) and the one- loop condition for 
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Figure 1: Phase diagram at A a = 0.06, A 2 = 0.002, A = -0.016 for a 8 3 
lattice. Error are smaller than the size of the marks. The points marked 
with triangles, for K\ = 0.33, are those for which we carried the detailed 
study of ISB. 



ISB eq. (|T3|) together imply 

A!>9^) 2 A 2 (14) 

and near the gaussian point K\ ~ k 2 ~ 1/3. Another condition was that 
Ai and A 2 had to be small in order to be near the gaussian point, but not 
too small because then it would be difficult to distinguish the various phase 
transitions one from the other. Moreover, it was desirable to have Ai as large 
as possible, in order for the stability condition eq.(^) and the ISB condition 
eq.([T3|) to be satisfied by a wide range of values for A's. The outcome of these 
considerations were the values Ai = 0.06 and A 2 = 0.002. The condition of 
stability for the potential then gives a lower bound on A of m —0.0219, while 
the one loop calculation predicts ISB for A < —0.012 ( assuming k\ ~ k 2 ). 
Inside the interval given by these two bounds, we have selected two values of 
A. The first is A = —0.016 which is well inside the stability region but not 
too close to the uncertain upper bound for the onset of ISB given by the one 
loop calculation. The second is A = —0.020: it is the one closest to the lower 
stability bound that we could use, without running into problems due to the 
nearby instability region. 
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In the simulations, we used a Metropolis algorithm combined with a Wolff 
single cluster method (the latter updates the sign of the fields), generally 
taking a ratio of 20 clusters every 3 Metropolis iterations (at the points 
where we found a first order transition the clustering was not efficient). The 
job was carried out on our RTNN computer, which holds 32 PentiumPro 
processors, for a total CPU time of approximately one month of the whole 
machine. The errors in the estimation of the observables have been calculated 
with the jackknife method. 

For each of the two triplets (Ai, A2, A), we started by drawing the approx- 
imate phase diagram in the plane (/cj, K 2 ) for T = 0, namely using symmetric 
lattices. This was done by simulating an 8 3 lattice and using hysteresis cy- 
cles to roughly identify the transition lines. Afterwards, in correspondence 
with each of these approximate critical points, we performed on the same 
lattice a simulation with much larger statistics obtaining a more accurate 
estimate of the transition as the maximum of the derivative along one of the 
axis of the Binder cumulant defined below. The resulting phase diagram for 
A = —0.016 is shown in fig. [I]: we observe that the k 2 ) plane is divided in 
four regions corresponding to the four possible cases of (fa) and (fa) being 
zero or different from zero. For A = 0, the critical lines would have been two 
straight lines each parallel to a coordinate axis. The interaction among the 
fields curved them in such a way that if one moves along, say, a vertical line 
parallel to the n 2 axis, starting from a point in the disordered phase, one 
encounters first the fa transition; if the value of k 2 is now further increased, 
one eventually crosses also the transition line for the fa field. This fact has an 
intuitive explanation: after the field 4> 2 has undergone the phase transition, 
it acquires a non-vanishing vev which, through the quartic negative coupling 
with the <pi field, acts on it as an effective negative mass term. Increasing k 2 
makes (fa) increase, until when the induced mass term for fa becomes large 
enough to cause a phase transition for the latter field too. A similar behavior 
is observed along lines parallel to the K\ axis. 

Having got an idea of the phase diagram, we passed to the accurate 
determination of some critical points. Since they had to be found with a 
high numerical precision, in order to clearly distinguish them from those of 
the asymmetric lattices, we could not afford to explore the entire (k\,k 2 ) 
plane, and thus we fixed once and for all the value of K\ = 0.33, and searched 
on that vertical line for the accurate critical values of k 2 , k 2 , corresponding to 
the transitions of <p 2 and fa. For A = —0.016 we simulated symmetric lattices 
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having N s = 8, 12, 16, 24, 32, while for A = -0.020 we had N s = 6, 8, 12, 16 
and could not use larger lattices due to the large correlation times. In any 
case, the results were very stable with N s and larger values were not needed. 
For each of the points and lattices that were simulated, the combination of the 
Metropolis algorithm with the Wolff single cluster method described above 
was repeated 1,250,000 times, taking measures every 5 iterations, while in 
the case of the first order transitions sometimes we made up to 5 millions 
iterations. 

The estimators used and the method followed to take the thermodynamic 
limit, in the case of second order phase transitions, were exactly the same as 



those used in ref. |22] and we refer to it the reader for the details. Here we just 
give a short review. To identify the critical points, we took the approximate 
values of k 2 given by the 8 3 lattice, and in correspondence with them we 
measured the Binder cumulants relative to that of the two fields which was 
undergoing the phase transition. We recall that the Binder cumulant of the 
field (bj is defined as: 



tt ^_ 3 ( M i) 
U Na;i {K 2 ) - - - , (15) 



where 



M i = \^ y E,M' ( 16 ) 



Having measured the Binder cumulants, we extrapolated them in a narrow 
k 2 interval around the simulation point using the spectral density method. 
The values of k 2 in the thermodynamic limit were obtained by looking at the 
intersections k 2 '(N s1 , N s2 ) among all possible pairs of curves U^ a>i {K 2 ) and 



using the following scaling law [28 



4(N S , bN s ) -4 = ^—N--^ , (17) 

where uo is the exponent for the corrections to scaling. On the symmetric 
lattices we used the exponents of the Ising model in three dimensions, namely 
v = 0.63, uj = 0.8, obtaining satisfactory fits. 

In the case A = —0.020, the shape of the phase diagram is very similar 
to that of A = —0.016, except for one very significant difference: while the 
transition of <p 2 is still of second order, that of (pi turned out to be of first 
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order (for A = —0.016 both transitions were of second order). To determine 
the precise location of this transition, we did not use the same method as for 
the second order ones. For first order phase transitions, general arguments of 
FSS do not apply and, to be specific, the crossings of the Binder cumulants 
do not coincide with the critical point. In this case, we searched for the value 
of k 2 , k>2(N s ), giving a maximum of the specific heat: 

where E\ is the Ki-energy of the field 4>x- 

^l = ^( E Ml,r + I *), (19) 
reZ 3 ,/i 

which is known to give a good signal for first order phase transitions. Using 
this estimator, the first order character of this transition appears to be very 
clear: fig. |] shows how the double peak structure, characteristic of first order 
transitions, becomes stronger as N s increases. The figure refers to symmetric 
lattices, but for all N t we have observed a similar behavior. This figure can 
be compared with fig. ^| corresponding to the <pi transition at A = —0.016 
and Nt = 3 as an example (here it is apparent the second order character). 
In fig. [| we show a typical Monte Carlo evolution for the first order case. 

Figure [5] gives an example of the fi^-energy and of the corresponding 
specific heat, extrapolated from the simulation point by means of the spectral 
density method. 

The critical coupling in the thermodynamical limit, i.e. k%, has been 
obtained from the measurements of K2{N S ) by means of a fit 

K2(N a ) - k% ~ N a - X ' v , (20) 

using for the critical exponent v its value for first order phase transitions 
in three dimensions, v = 1/3. The fits turned out to be good in all cases, 
supporting the first order character of the transition. 

Having determined the phase diagrams at T = 0, we passed to the study of 
the finite temperature case, by simulating asymmetric lattices with N t < N s . 
The values of N t that we considered have been N t = 2,3,4 for both values 
of A and Nt — 5 in one case. The values of N s used were the same as for 
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(18) 




Figure 2: Hystogram of E 1 for A = —0.020 on the X transition (the dotted, 
dashed and solid lines refer to the 6 3 , 8 3 and the 12 3 lattices respectively). 




Figure 3: Hystogram of E 1 for A = —0.016 on the X transition (N t = 3, dots 
for N s = 12, dashes N s = 16, daashes iV s = 24, solid iV s = 32). 
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Figure 4: Monte Carlo evolution of Ex for A = —0.020 on the <f>i transition 
for the 12 3 lattice. 

the T = case, except that for A = —0.020 we have sometimes gone up to 
N s = 24. The order of the transitions turned out to be always the same as 
in the corresponding T = case; in particular, for all N t 's we have found 
the first-order transitions described above for the lower value of A. The 
techniques used to take the thermodynamic limit N s were the same used in 
the T = case, except for two differences. Since the scaling parameter was 
N s , while N t was fixed, according to the hypothesis of dimensional reduction 
and universality, we used the exponents of the Ising model in two dimensions, 
namely v — 1 and to = 4/3. Similarly, in the case of the first order phase 
transition encountered for A = —0.020, in eq. (|20|) , we used the value of v for 
two dimensions, namely v = 1/2. In both cases we got good fits. 

As far as ISB is concerned, the most important result of the simulations is 
apparent from figs. [] and 0: for a constant K\, they show that, starting from 
the symmetric lattices (the points with 1/N t = 0), k\ for the transition of 
the (f>2 field increases monotonically when Nt is decreased for both values of 
A considered; this means that the critical points for Nt finite shift deeper and 
deeper in the ordered region of the T = model, i.e. raising the temperature 
disorders the system, as it happens in normal cases (see for a comparison 
p2]j ) and contrary to what is required for ISB to occur. Though we have 
considered a single value of K\ in our simulations for finite NtS, we do not 
expect a different behavior at other points. If for some other value of k% the 
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Figure 5: E 1 and C\ for A = —0.020 on the 0i transition for the 12 3 lattice 
by using the spectral density method. The full mark denotes the simulation 
point, the empty ones are extrapolations thereof. 
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Figure 6: k 2 c as a function of 1/N t for (Ai = 0.06, A 2 = 0.002, ; 
—0.016, Ki = 0.33). The points with 1/N t =0 refer to the symmetric 
tices. 
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Figure 7: k 2 c as a function of 1/N t for (Ai = 0.06, A 2 = 0.002, A = 
—0.02, K\ = 0.33). The points with 1/N t =0 refer to the symmetric lattices. 
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critical A^-line passed below that for the T = theory, they would have to 
cross at some point, which would be a very strange fact. 

With such a result, there is no need of looking at the continuum limit 
to rule out ISB. Recalling the comments at the beginning of this section, 
it is obvious that the CCP's associated with continuum theories having a 
symmetric vacuum at T = cannot cross in any way the 2 critical lines 
of the finite N t lattices, for the latter lie in the ordered region of the phase 
diagram for the symmetric lattices, and not in the disordered one. Then this 
study indicates that Inverse Symmetry Breaking does not work in this model. 
However, it is not conclusive about Non Symmetry Restoration because we 
do not know if the temperature for which the symmetry is restored diverges 
or not in the continuum limit. 

5 Conclusions 

Inverse Symmetry Breaking and Symmetry Non-Restoration have recently 
attracted much interest, mostly due do their potential implications for cos- 
mological scenarios. They are based on the possibility that in multiscalar 
models some Debye masses can be negative at high temperatures, for suit- 
able choices of the couplings. 

In this paper we have presented the results of a Monte Carlo study of ISB 
in a simple model for two real scalars interacting via a \4>l4>l coupling, in 
three euclidean dimensions. The reasons for carrying such a study is that in 
three dimensions perturbations theory cannot provide any information con- 
cerning these phenomena, as it is discussed in Sec. 3. In order to have a fully 
non-perturbative approach to these phenomena, we have thus simulated this 
scalar model on the lattice and found no trace of ISB. It turned out that for 
all values of A that we considered the size of the disordered region of the phase 
diagram increases, when the extension N t of the lattice in the euclidean-time 
direction is decreased, as it occurs in normal systems, instead of decreasing 
as required for ISB to occur. Heuristic arguments tell us that this conclusion 
should be true for the whole parameter space. Another interesting result of 
our simulations is that, for large negative values of A, one of the transitions 
becomes of first order. 

It would be desirable to know if other non-perturbative analytical meth- 
ods can reproduce our findings. Among the possible approaches we mention 
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here the one based on the Average Effective Action f26|, or the CJT for- 



malism [29], that have already been used to study the thermal behavior of 



two-scalar models in four dimensions |T(| [T7|, |23J . 

We are currently simulating the model considered in this paper in four 
dimensions, which is physically the most relevant case. Our efforts are being 
focused on two issues, both very relevant for Cosmology: one is the study of 
ISB and SNR and the other is the existence of first order phase transitions 
at finite temperature. We hope to report soon on this work. 
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